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The symmetries of the DNA double helix require a new term in its linear response 
to stress: the coupling between twist and stretch. Recent experiments with torsionally- 
constrained single molecules give the first direct measurement of this new material pa- 
rameter. We extract its value from a recent experiment. We also present a very simple 
microscopic theory predicting a value comparable to the one observed. Finally we sketch 
the effect of constrained twist on entropic elasticity of DNA arising from the connection 
between Link, Twist, and Writhe. 
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1. Introduction 

The idea of studying the response of DNA to mechanical stress is as old as the discovery 
of the double helix structure itself. While many elements of DNA function require detailed 
understanding of specific chemical bonds (for example the binding of small ligands), still 
others are quite nonspecific and reflect overall mechanical properties. Moreover, since the 
helix repeat distance of £q ^ 3.4 nm involves dozens of atoms, it is reasonable to hope that 
this length-scale regime would be long enough so that the cooperative response of many 
atoms would justify the use of a continuum, classical theory, yet short enough that the 
spatial structure of DNA matters. Since moreover various important biological processes 
involve length scales comparable to £q (notably the winding of DNA onto histones), the 
details of this elasticity theory are important for DNA function. 

Recently, techniques of micromanipulation via optical tweezers and magnetic beads 
have yielded reliable numerical values for the bend stiffness from the phenomenon of 
thermally- induced entropic elasticity (Smith et al., 1992 Bustamante et al., 1994; Volo- 
godskii, 1994; Marko et al., 1995), as well as the direct measurement of another elastic 
constant, the stretch modulus, by exploring the force range 10-50pN (Cluzel et al., 1996; 
Smith et al., 1996; Wang et al., 1997). Significantly, the relation between bending stiffness, 
stretch modulus, and the diameter of DNA turned out to be roughly as predicted from 
the classical theory of beam elasticity (Smith et al., 1996), supporting the expectations 
mentioned above. 

Still missing, however, has been any direct physical measurement of the elastic con- 
stants reflecting the chiral {i.e. helical) character of DNA. Recent experiments with tor- 
sionally constrained DNA have permitted the determination of one such constant, the 
coupling between twist and stretch (Strick et al., 1995; Marko, 1997; Kamien et al., 1997). 
This coupling may be relevant for the binding of the protein RecA to DNA, which stretches 
and untwists the DNA (Stasiak et al., 1982). We will explain why this term is needed, 
extract its value from the experiment, and compare it to a the prediction of a simple 
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microscopic model to see that its magnitude is in line with the expectations of classical 
elasticity theory. Finally we will briefly sketch how to understand another phenomenon 
visible in the data, the effect of constrained link on entropic elasticity (J.D. Moroz and P. 
Nelson, in preparation). 

2. Experiment 

DNA differs from simpler polymers in that it can resist twisting, but it is not easy to 
measure this effect directly due to the difficulty of applying external torques to a single 
molecule. The flrst single-molecule stretching experiments constrained only the locations 
of the two ends of the DNA strand. The unique feature of the experiment of Strick et al. 
was the added ability to constrain the orientation of each end of the molecule. 

We will study Fig. 3 of (Strick et al., 1995). In this experiment, a constant force of 
8pN was applied to the molecule and the end-to-end length z^^^ monitored as the terminal 
end was rotated through ALk turns from its relaxed state (which has Lkp turns). In this 
way the helix could be over- or undertwisted by as much as ±10%. Over this range of 
imposed linkage z^^^ was found to be a linear function of a: 

£ = const. — 0.15a where a = ALk/Lkg and e = (^^^^/z^q^. q) — 1 . (2.1) 

Thus a is the fractional excess link and e is the extension relative to the relaxed state. 
Eqn. (p.lj ) is the experimentally observed twist-stretch coupling. 

The existence of a linear term in ( |2.1|) is direct evidence of the chiral character of 
the molecule, and its sign is as expected on geometrical grounds: untwisting the molecule 
tends to lengthen it. Still geometry alone cannot explain this result. Consider the outer 
sugar-phosphate backbones of the DNA. Suppose that the twist-stretch phenomenon were 
due to the straightening of these helical backbones while they maintained constant length, 
0.6 nm per phosphate, and constant distance 0.9 nm from the center of the molecule. Then 
since each basepair step is h = 0.34 nm high, the circumferential length per step is = 
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-\/.62 — .342 nm. The corresponding twist angle per step is given by ^ = (^^Z^) /-Qnni = 32°, 
roughly as observed. Supposing now an extension by A/i//i = e, we find an untwisting 
by a = 56/6 = const. — e/2.0, quite different from what is observed, eqn. ( p.l| ). We must 
seek an explanation of the experimental result not in terms of a geometrical ball-and-stick 
model but in the context of an elastic response theory. 



3. Simple Model 

We will begin by neglecting bend fluctuations (see below). A straight rod under 
tension and torque will stretch and twist. We can describe it by the reduced elastic free 
energy 

= ^ [Ca' + Be' + 2Dea] - fe . (3.1) 

Here C is the twist persistence length, B ^ 1100 pN/c^Q /c^T ^ 78 nm is the stretch modulus 
(Wang et al., 1997), and D is the desired twist-stretch coupling. L is the relaxed total 
length, Wq = 271 / Iq = 1.85/nm, and the reduced force / = 8pN//c^T ^ 1.95/nm in the 
experiment under consideration. For a circular beam made of isotropic material the cross- 
term D is absent, since twisting is odd under spatial inversion while stretching is even. For 
a helical beam, however, we must expect to find this term. 

We now minimize F with respect to e at fixed force with an imposed constraint on 
the overtwist a to find 

e = e^^o - iD/B)a . (3.2) 
Comparing to ( |2.1j ), we obtain the desired result: D = 12 nm. 



4. Bend Fluctuations 

We have discussed the term linear in overtwist a in ( |2.1|) . For the highest-force curve 
at 8 pN this is the dominant effect. At lower forces, however, it is quickly overwhelmed by 
an effect symmetric under a —cj, which we have so far neglected. This effect is due to 
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the coupling between applied overtwist and thermal bend fluctuations. We now sketch a 
simple, though imprecise, analysis of this effect. The full analysis is qualitatively similar 
(J.D. Moroz and P. Nelson, in preparation). 

Since in this section we want to study nonchiral, low-force effects, we will revert to 
a model of a fixed-length cylindrical rod with bend and twist elasticity. We will consider 
small deviations from the unstressed state of the rod, which we take to run along the z 
axis. Initially we paint a straight stripe on the outside of the unstressed rod. To describe 
the deformed rod, we find at each point a triad of unit vectors {Ej(s)}, where E;^ is the 
tangent to the curve determined by the rod centerline, £3 -L is the normal vector from 
the centerline to the stripe. Eg = Ej^ x Eg, and s is arclength. Let 77 = o^qCt be the 
imposed excess helix density, and define the convenient reference frame 

^^{s) = xcos(?7s) + y sin(?7s) ; ^2^^) = — xsin(?7s) + y cos(r7s) ; 63 = 2. 

We can now describe the deformed rod by three small variables: the projection tj^ = 
'^i^i + ^2^2 of -^3 plcine and the angle ip between x and the projection of Ej^ to 

the xy plane. We propose to expand the elastic energy to quadratic order in these and 
thus find the thermal fluctuations in harmonic approximation. In terms of tj^, <^ we find 

El = (1 - - |^')ei + - (ti + ^2^)2 

E2 = -((^ + tit2)ei + (1 - 1*2^ - \V^)^2 + i-h + 

E3 = (1 - \ti - W> + ^lei + ^2^2 • 
We may now differentiate with respect to arc-length to get the body-fixed angular velocities 

^1 = E3 ■ E2 = -7-2 + Wl, ^2 = El ■ E3 = Ti + <pT2, O3 = E2 ■ El = ?7 + - |?7(ti2 + 

^2^) +^1^2' where we abbreviated = i^— 77^2, ^2 = ^2 Note that the formula for 

the Twist, fl^i J^^^ ^ simple derivation of Fuller's formula (Fuller, 1978) for the Writhe 
of a nearly-straight curve. 

Our formulas become very compact if we introduce the complex variable T = (x -|- 
iy) • t^. We then have + = I'J'P and r^^ + T2^ = |Tp. Finally we expand in 

4 



Fourier modes: T = J2q tt^e"^'^ and similarly with Substituting into the elastic energy 
E/k-^T = ^ J ds [A{fli^ + 02^) + C'^s^] yields the harmonic elastic energy 

E/k^T=lJ2[^^'-C^^ + f] I^Z + ^C'E^'I^/- (4.1) 

We have introduced the applied stretching force / = f /k^T. 

The physics of this formula is clear: twist fluctuations decouple from bend fluctu- 
ations, but the imposition of nonzero net overtwist 77 creates a new crossterm in Og^, 
first-order in s-derivatives, and this crossterm affects the long-scale bend fiuctuations. In- 
deed completing the square in (|4.1[) shows that the effect of 77 is to reduce the effective 



tension from / to / — ^fj-{Crf)'^ . This effective reduction is what makes the relative ex- 
tension plummet as the overtwist rj = Wq cr is increased at fixed /. This effect, combined 
with the intrinsic twist-stretch effect from the previous section, explains qualitatively all 
the phenomena in the region of the experiment where linear elasticity is valid. A more 
precise version of this calculation also affords a direct determination of the value of the 
twist stiffness C (J.D. Moroz and P. Nelson, in preparation; C. Bouchiat and M. Mezard, 
in preparation.). The high-force regime studied here is free from some of the difficulties 
of the Monte Carlo approach (Vologodskii et al., 1979; J.F Marko and A.V. Vologodskii, 
submitted); in particular, there is no need for any artificial short- length cutoff. 



5. Microscopic Model 

The elastic theory in §3 was very general, but it gave no indication of the expected 
magnitudes of the various couplings. To gain further confidence in our result, we have 
estimated the expected twist-stretch coupling based on the measured values of the other 
elastic constants and geometrical information about DNA (Kamien, 1996; R.D. Kamien, 
T.C. Lubensky, P. Nelson, and C.S. O'Hern, submitted). We used a simple, intuitive 
microscopic picture of DNA as a helical rod to show how twist-stretch coupling can arise 
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and get its general scaling with the geometric parameters. The model shows that the value 
of D calculated above is reasonable. 

6. Conclusion 

We have pointed out a strong twist-stretch coupling in torsionally-constrained DNA 
stretching experiments, evaluated it, argued that it reflects intrinsic elasticity of the DNA 
duplex, and shown that the value we obtained is consistent with elementary considerations 
from classical elasticity theory. We also showed how the interplay between twist and 
writhe communicates a constraint on link into the entropic elasticity of DNA, as seen in 
experiment. 

We would like to thank D. Bensimon, S. Block, and J. Marko for their help and for 
communicating their results to us prior to publication, and W. Olson for discussions. RK, 
TL, and CO were supported in part by NSF grant DMR96-32598. PN was supported in 
part by NSF grant DMR95-07366. 
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